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Abstract 

We explore further the classical QGP using methods from classical liquids. The partition function 
of an ensemble of SU{Nc) colored charge spheres is constructed. We analyze it using a cumulant 
expansion (low density) and a loop expansion (high temperature) after resumming the Debye screen- 
ing effects. The pertinent free energies are derived in both limits and compared to recent molecular 
dynamics results. 
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I. INTRODUCTION 



Plasmas are statistical classical/quantum systems involving charge constituents. Notable 
are electromagnetic plasmas whereby the underlying constituents interact through long range 
Coulomb fields. The simplest theory of electromagnetic plasmas is the One Component 
Plasma (OCP) whereby the constituents are like (negative) charges embedded in a uniform 
and neutralizing unlike (positive) background. A number of analytical approaches to the 
OCP exist in the form of diagrammatic or field theoretical methods Q]- These formal 
approaches form a useful theoretical corpus for understanding ionic liquids. 

Ionic liquids are characterized by a short range repulsive core in addition to the long 
range Coulomb potential. The core is a scalar with a range of the order of the interparticle 
distance, while the Coulomb attraction is vectorial with infinite range. Reference liquids 
are described by a repulsive pair potential w{x) = oo for x < a and w{x) = for x > a 
where a is the diameter of a hard sphere. The Coulomb potential is generally viewed as a 
perturbation added to the repulsive core. Many models of ionic liquids have been developed 
to accommodate arbitrary charges and cores such as the Primitive Model (PM). Spinoffs are 
the Restricted Primitive Model (RPM) with opposite charge pairs and the same core, and 
the Special Primitive Model (SPM) with arbitrary charges and the same core. 

The Classical Quark Gluon Plasma (cQGP) as developed by Gelman, Shuryak and Zahed 
can be regarded as an example of an SPM model in ionic liquids, albeith with non-Abelian 
colored charges [s]. The nature of the core in the CQGP is quantum mechanical and thus 
assumed. Detailed molecular dynamics simulations of the cQGP sl have shown a strongly 
coupled plasma for T = V/K ^ 1, i.e. whenever the potential energy is of the order of the 
kinetic energy. The classical and colored cQGP maybe in a liquid state at moderate values 
of r prompting us to use methods for classical liquids to analyze it. 

This paper is the first of a series of sequels to jsl] to develop analytical methods to address 
the many-body dynamics of the cQGP. In section 2, we define the partition function of the 
cQGP as a classical but colored liquid. In section 3 and 4, we work out the partition function 
in the low density limit through a cumulant expansion after resumming the Debye screening 
effects. In section 5, we discuss a high temperature expansion. In section 6, we detail a 
diagrammatic or loop expansion for the cQGP that is justified at high temperatures. We 
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carry the expansion to three loops. In section 7, we unwind the free energies both in the 
loop and density expansions. In section 8, we compare the excess free energies to recent 
molecular dynamics simulations. Our conclusions are in section 9. Useful definitions and 
color integrations are given in the Appendices. 

II. THE GRAND CANONICAL PARTITION FUNCTION FOR CQGP 

The cQGP has been defined in jsj. It consists of classical particles of like-mass m with 
3-position Xi{t) and 3- momentum Pi{t) with A^"^ — 1 adjoint colored charges Qf{t) and fixed 
Nc — 1 Casimirs. Because of the constraints, the color variables are Darboux's type as 
summarized in the Appendix, (see also ^\). Particle motion in the cQGP is treated classically 
through 

= gQ'^Frix.^pUt) 

= -9r''p^{t)Al{x.)Qm (II.l) 

where we have trivially covariantized the notations. The repulsive effects of the core is 
subsumed jsl. A'^ is the gauge field on the ith particle due to all other particles, F^" is 
its field strength, and f"^^ the pertinent structure constant for SU{N^. The last relation 
is Wong's equation [5]. The magnetic contribution to the Lorentz force is suppressed by 
v/c and dropped in the electric cQGP. cQGP simulations with both electric and magnetic 
charges can be found in P]. 

The Hamiltonian for the cQGP reads 

Inia 2 47r -7^, ■ \ Xai — X'l A 

where the core potential is now explicit. The double sum is over species a (to be set to 3 
below) and particle index i = 1, .., A^. In terms of the color charge density 



P''{^=Y.QaA^-^a,) (II.3) 



3 



the grand partition function for the cQGP is 



^ = E / n Tn X{Wa4'ra,ina) R ^xp " / d?rd\' p^^ {r)v{r - i')p''{r') (II.4) 
with 

f(r — 7^) = — --^ — — (II. 5) 

^ ' 4:Tx\r-r' \ ^ ' 

and the free particle density 

ua = 9a^e^^- . (II.6) 

Qa is the degeneracy factor and /i" the chemical potential of species a. The factor of A^'^ 
is the thermal wavelength resulting from the momentum integral over phase space. The 
repulsive core potential is set to be 

^ \ oo (\ f — r' \< a) 
w{f-f^) = { ^' ' ' (II.7) 

[O (|f-F|>or) 

with a core size a to prevent classical collapse. In reality, the core emerges from Coulomb 
repulsion between like-particles or the quantum uncertainty repulsion between unlike/like 
particles. With this in mind, we can rewrite (1II.4P in the form 



^ 1 Na N^-l 



Na a ^"a- a,i 

X exp ^ — cP'rd^r' p"{f)v{f — f) p°' (r')^ e^^° exp ^ — | / d^rd^r' p(r)w{r — r')p{r' 

(II.8) 

with the number density p(r) = J2a,i S{f— ra,i)- Wq = w(0) is identified with the self-energy. 

The repulsive core (111.71) acts like the potential between two hard spheres with diameter 
R = a. Thus, we may treat each particle classically as a rigid sphere of diameter R with 
uniform color charge on the surface. 



Qq{S) = -%-^5{\x\-R/2) (II.9) 
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The repulsive core flll.7p is now generated classically as the Coulomb repulsion between the 



rigid spheres 



7|. For two spheres located at r and this is 



W{r — r)=[ dxdyq{\ r — x |)— f — y \) . (11.10) 

J \ X — V \ 



y 

The grand canonical partition function follows in the form 



r 1 



-({N^-l)Wo^-VHS 



Na « a,i 

X exp ( - f g e5"^ / d^rd^r'p''{r)W{f- F)p"(r^)) (11.11) 



where the Coulomb self energy is now Wo = W{0)/2 and 7a 



III. LIQUID OF COLORED HARD SPHERES 

To analyze (III. lip we first consider a dilute cQGP with Nc colors. Performing (A*"^ — 1) 
times the Sine-Gordon transform on (III. lip yields 



exp ( - f-^ I dfd^p''{r)W{r- r )p"(r )) = (exp ( ^il^^f'^ I dfp''{r)(j)''{r)))w 

2 47? ex. <y. 47r J 

(111.1) 

where the averaging on the RHS (right hand side) is carried using the measure 



/[n? d(p%- ■ ■) exp ( - I '"^ / drdf(i)''{f)W-\r~f)(j)''{f, 
To proceed further, we define 



= I3p, = + {Nl - 1)71^0 na = nae(^'-i)^^»e'(^-)"S"'"S-''^-'^"(^^-^) (III.3) 

as the renormalized chemical potential and species densities respectively. In terms of which 
the grand canonical partition function is now 
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r 1 

^ = E / n TFT X{{dQa,d?raAna)w)e-"'"' = {Zhs)w (111.4) 

Q-'"Hs ig short for the Gaussian measure in (1111.2^ . (1111.2^ can be Taylor expanded around 
the mean density (see below). The result is 



\ZhsJ t^inl a,,^..,aj 



5™ InZ 



ni 



'^^ai(^ai) ■ ■ ■^^a„(^a„) nai(raj=na i=l 

If we were to assume that the species chemical potentials are all the same fiai 
Ha, then in the free particle case 



(III.5) 



A^o = /ia + ln(3), no = Sn^ = 3A 



(III.6) 



We have set the particle species to 3 to account classically for quarks, anti-quarks and gluons 
all of equal constituent thermal mass m for simplicity. In terms of (1111.60 the integrand 
in flIII.3 ) can be identified with the classical correlation function /ig"^ of a liquid of hard 
spheres [7| 



n. 



with the number density 



n 

|r4"^(ra„ ■ ■ -r-J = p-X\ra,r ■ -r.J (III.7) 



Po 



1 dlnZns 



V duo 

Thus, the partition function for a colored liquid of hard spheres now read 



(III.8) 



In 



{exp{-U[r]))w = Afw' / [Y[dr]exp{-S[r]) 



with 



(III.9) 



Mw= I [\{d(t>^] exp ( - 1 / dTdf'ct)'^{T)W-' 



r — r ]0 [r 



S[<P"] = ^ E / dfd7^(l)'^{f)W-^{f- r')0"(r') + ^ f/^^-^^^ 

^ a n=l 

^- ai,a2,--,an i=l 

IV. CUMULANT EXPANSION 

To perform a cumulant expansion at low density, we add and subtract the Hartree part 
in ?7[0] through 



f/n = (f/r] - t/ori) + f/o[0i (iv.i) 



with 



^PoET"E/*>°(0% (IV.2) 



a a 



yields the identity 



(exp(-f/r]))H. = :^(exp(-f/r]))x (IV.3) 
with =W~^ + po Sa7al- The normalizations in (1IV.2I) is explicitly 

^ = exp ( - ^{Nl ln(l + Po E 7aW^(^))) (IV.4) 

The Fourier transform of the core W{q) is related to the Fourier transform X{q) as 

with 



All in all the partition function for a colored liquid of hard spheres with their color charges 
smeared on the surface is 
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(IV.7) 

with 7i[0"] = f/[0"] — f^o[0'^] and (■ ■ ■)x,c denotes a cumulant average. At low density, the 
grand partition function is dominated by the first cumulants 

-i^ = -i^ + ™, + -. + 0(p^) (IV.8) 

with 



"^•i = - 1) / ^ + W E7<.W(?)) 

H = f/i [01 + U2 [01 - f/o [r ] (IV.9) 



The leading or zeroth order cumulant is known from hard sphere liquid analysis 



= -p; - f - ^-v?<^'' + o(p?) (IV. 10) 



with pg = 3po- The first cumulant reads 



w,-{Nl-l)p,Y.^aW{^) 

a 

= -{Nl - l)\V^PonY.laf" + {Nl - l)^vrpo^(E7a)^i? 

-(iVf - l)\ir"'p.nY.laf"R' + O(P^) . (IV.ll) 

a 

The Gaussian averaging in the second cumulant can be carried out. The result is 

^2 = -{Nl - l)^E7a - \{Nl - l)poE A.0 + Po(3 - Ee-^(^-^)^-) 
-f ( E 7a(l - e-^(^-i)^-)) ' / rffX^(r) - ip^(3 - E e-^(^-^)^-)'/.o(0) 
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a,a' 

x(e-^l^^«E.QSQ:, _ 1 _ l(Ar^2 _ i)^„^„,X^(r))(/if^(r) + 1)) (IV.12) 
with = ^aX(0) — 2^a/R ^nd ho(0) the 3D Fourier transform of hQ'\r), 

A„o = 7a^(0) - 2^ = 7„( - «o + - + 0{4)^ (IV.13) 

and 

7„X(r) = 2^ - - 7„«:o + (r < a) 

Collecting all the results up to the second cumulant in powers of the original activity 

'^o(= A) we have for Nc = 2 



- ^ = - 2V^^'''(E^«)^/' - A^(^-E^aE^a + ^2') + 0(A^/^) (IV.15) 
up to order A^/^, with 



S2' = ^ E (54e' /^'' ^c^t + e'^ (-2 - 2e - 3e^ - 9e=^) + e-^^(2 - 2e + 36^ -9e')) (IV. 16) 

a, a' 

where = e = y^e^e^ and A = cr^A. We have set h'^\r) = —1 for r < cr and h^\r) = 
otherwise, and used 



The correction of order A^/^ is straightforward but tedious. For any Nc we obtain 
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a a ^ a a 

HN! - l)(E7a)^(E7a)^27ria^ - ^iN^ - l)^(E7D(E7a)^2vria 



3 

,712 



-v^(E7a)^/^Er^^ 

(7a + 7a'o) . 



1. 1 



-Tsinh(7,,,(iV,' - 1)-) - rcosh(7„,,(iV,' - 1)-) 

1 



+ ^':- T' sinh(7„.KiV,^ - 1)-) + -iNl - 1)(7„ + 7a') + lala'iK - 1)- 
^laa'- r 2 r 

-\{la + la^lala\K - 1)^1) (IV.18) 

The left out integration is divergent. For A^c = 2, the result is logarithmically divergent 



^vri($:7.)nE7D^ln(-) (IV.19) 

where A is an infrared cutoff in the radial direction. This result is in total agreement with 
the result obtained earlier in Sj using a method that did not account for hard spheres and 
the smearing of the color charge. For arbitrary A^^^ the order A^^^ correction is 



-(iVe^-lf^(E7a)nE7D 

a a ^ ^ 

HN! - l)(E7a)^(E7a)^2vrfa^ - W - l)^(E7D(E7a)i2vrfa 

a a ^ a a 

+ ^-i(E7a)nE7D^ln(^) 



(iVf-lf^E7.^(E7a)^ 



(IV.20) 



V. HIGH TEMPERATURE EXPANSION 



Another useful way to analyze (jIII.9tiIII.10p is through a high temperature expansion 



which parallels the cumulant (yirial) expansion. For that we order the last part of f/„ [0°^] in 
(Illl.lOj) in powers of /5 as in S]- For that we define 7' = Pg'^/Air so that 



1 = 7a(iV,2 _ ^ Z7'^ 5:Q:0"(r) - 7'^ E QX'ririr'if^ 



R 
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+ E QaQa'QfQf0"(r)0"'(r)0"''(r-)0"'''(f) + O(/?t) (V.l) 



oi,a'a",a"' 



The rules for color averaging are summarized in Appendix B. For instance, using flB.7p and 
(IRTTD we have 

hUimx = \poY.^ao{Nl - 1) - \poY.^lM - 1)' + 0{^') (V.2) 
in leading order. In particular, the second cumulant at high temperature is 



^ a ^ a 

-^Po(E Aao)^(iV.^ - 1)^ - -ANl - l)p,\Y.laf j d'rhf\r)X\r) 

+ \jf^'i''iCD' j d'r{hf{r) + l)X^(f) + . (V.3) 

The last term contributes only for > 2 with ^3 being the cubic Casimir say for SU(3) as 
detailed in Appendix B. 

After collecting all terms and replacing X{r) with ^e~'^°^ in the limit i? — 0, we finally 
obtain 



InZ 



V V 



^-^(iV.^-l)pf(E7.P 

- l)V^(E7a)^ / d'rh^^\r) - 2n{N^ - l)pg(E7a)^ / d'rh^^\r) 
(iVf-l)V^(E7D(E7a) 



TT 

2 

Below we show that this result yields a high temperature free energy for the cQGP that 
is identical to the one following from the loop expansion with an infinite core ( 111.7^ with 
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VI. LOOP EXPANSION 



The loop expansion of (III.SP is best captured by reorganizing the expansion around the 
Debye-screened solution. This expansion is identical with the high temperature expansion 
of the hard sphere liquid in the limit of zero size spheres. The finite size case will be derived 
by inspection. With this in mind, we can perform and perform the Hubbard-Stratonovitch 
transform on the colored interaction part of (111.81) . 

exp (^"^f^ J drdr* p°'{r)v{r — r)p"{f^)j 

= (^/3 det{v~^)^ ' J exp - ^ y c/f0"(r)t;~V"(O) exp J rffp"(f)0"(f)^l.l) 
and similarly for the core part in (111.81) . The partition function for the cQGP reads 

Z = i^pdetiv'-^)j i^detiw-^)j j n (VI.2) 

with the induced action 



S = J2 - [ c^r0"(f)t;-V"(^) + - / t^f^(r)M;-V(^) 

- E / ciQaC^rn,e*'^S^"«"<^"(^^+*^(^^)+^-° (VI.3) 

a •' 

Here wq is the divergent self-energy. For simplicity the colored particles are point-like 
throughout. 

If we introduce the screened Coulomb potential 

VDH{r) = P-^e-^^^ (VI.4) 
with the squared Debye wave number, 

c ^ a a 

then the induced action (1VI.3P can be split into a screened part 5*0 (quadratic in the fields) 
and an interaction part Si (rest). Specifically, 
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N?-l ... n ,.2 

" a, a' 

= 5o + 5, (VI.6) 

where we used the color normahzation / dQ = 1 and the color averaged squared Debye wave 
number, 



,2 . N'^-l 



kI^ I dQX = -I— / dQ,Y.^a/3 E 

T I dQaY. ^aP{N! - l)C2a = Q^Y. ^af^^^a (VI. 7) 



iV2 

C u. 

Here C2 is the quadratic Casimir (^2) divided by (A*"^ — 1), 

j dQaQ:Q:' = C2.5""' = ^^1^5""' (V1.8) 

A. One-Loop 

The screened Debye-Huckel partition function follows by setting S*/ = in the induced 
action. The corresponding partition function is then 

Zo = exp (vYna) [det (1 + -^,^1)) (VI.9) 
The argument of the determinant is the inverse screened Green's function 



^{ -V^ + Kl)G{f-f') = S{f-f^) (VI. 10) 



which is 



Gif- r') = / / (VI.ll) 

The apparent singularity for coincidental arguments can be handled by dimensional regular- 
ization [lo| . 
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lim ^^(O) = ^(0) = -Ka^ (VI.12) 

with n being the spatial dimension. The determinant in flVI.91) can be calculated by standard 
methods. The identity 51ndetX = Tr{X~^5X) yields 

Indet (1 + -X-kI) = -GM'^V (VI.13) 

so that 

Zo = exp (viY^ria - — ^^^(O)) (VI.14) 



Using ( 1VI.7[) and ( ]VI.12[) . we obtain the screened one- loop result as 



Zo = exp [Vj2nail + {N^ - l)g^f3C2a 



127r' 



exp ( V ^ na{l + g'^(3q2a 



12^' 



(VI.15) 



B. Two-Loop 



Higher order loop corrections follow from 



Z = ZoU exp / dndr.^Gin - r.) 



exp (-5,(0")) 



(VI.16) 



00=0 



.2(3 J 5<p^ ' '6<p° 

by insering higher G's through Sj as recently discussed in [10] for the Abelian case. The non- 
Abelian is noteworthy in many respects as we note below. A typical two-loop contribution 
is shown in Fig. [H Its contribution is 



FIG. 1: Two loop contribution from one particle (IVI.lTp 



^^P\Nl - lfY.^aC2a'j dfG\0) (VL17) 
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The color factors follow from the identity (1B.7P discussed in Appendix B, i.e. 



n n n 

(VI.18) 

with A = 0, B = l/15gi for SU(2) and A + 3B = 3/80g| for SU(3). Although SU(2) and 
SU(3) involve considerably different integration measures, the overall result for Fig. [1] is the 
same. Additional contributions to the colored partition function are shown in Fig. [2|, which 
contribute 

"a I \ "a 



FIG. 2: Two loop contribution from two particles (|VI.19p 



1 



-qs J dfdr'G^ - ^) XI 

o. a' 

^^JdfdfG{f~f')G\0)Y.naY.na'^ E d^^^^d^^^ (VI.19) 

a' «,/3,T 



2-3 

(3' 



2 ^ . 

a' Q!,/3,7 

Here is the cubic Casimir following from the identity 



5: d^p^Q^Q^^Ql = 53 = ^ E d^p,d^^'< = J dQ^qs = E / dQado^MQ^aQZ 

a,(3,^ a,[3,-y ct,/3,7 

/ dQaQMQ: = \d''^^ (VI.20) 



as detailed in Appendix B. The cubic Casimir exists only for Nc > 2 and vanishes identically 
for SU(2). This is clear from the fact that the contributions in Fig. [2] involve 3 colored 
vertices. Also note that only the irreducible graph of Fig. [2] contributtes to fIVI.lQP since the 
reducible graph averages to zero by color integration through the identity 
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(^""/3^/377 = (VI.21) 

oi,l3fr- 

Unlike the Abelian case discussed in [lOj, where tadpoles and disconnected contributions 
abund, the non-Abelian case has none of these thanks to the color integrations. Also, the 
effects of the hard core at two loop can easily be recalled by noting that in the 2-particle 
channel under consideration the distance of minimum encounter is a. So the radial integra- 
tions should be limited to a < r < cxd to account for the hard core. With this in mind, the 
2-loop contribution to the partition function reads 



In Z In 1 



where the last term is only present for SU(3). Using flVl.151) and flV1.12p we get 



(VI.22) 



^ a '-' a 

a a 



(V1.23) 



2-3! Mvr' 4 

Here -Ei(x) is the exponential integral(— x) = -E'i(x) = ^dt), which is logarithmi- 



cally divergent at short distance. It is made finite by the hard core potential. 



C. Three Loop 

A partial three loop analysis will be carried out in this section. There are in total 8 
diagrammatic contributions at three loop that can be organized in terms of the particle 
density: 1 (one-particle density); 3 (two-particle density); 4 (three- and four-particle density). 
They will be considered sequentially. 

The three loop contribution stemming from the one-particle density is shown in Fi. [31 It 

is 

- dfG%0)q2a' (VI.24) 
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FIG. 3: Three loop contribution from one particle (|VI.24|) 

The vertex involves 6 color charges which are integrated with the help of the identity (see 
Appendix B) 

+d"^^d^"^ + d'^^'d'^^^ + d'^^^d'^^' + d'^^'d'^^^ + d'^^^d'^^' + d^'^df^^^) 
+b(5''^5^^5'^ + + r'^^^^^^^ + r^^/^^fj^^ + 5"^ 5^' 5^^ 

^^ae^p&pC + + d'^^d^^S^' + + (VI.25) 

with A = 0,3 = -^ql for SU(2) and A = + f 5 = f - §ql for SU(3). Note 

that this contribution is similar to one-loop for both SU(2) and SU(3) despite the differences 
in the contributions and the color averaging. 




FIG. 4: Three loop contribution from two particles (1VI.26P 

The three loop contribution stemming from the two-particle density are shown in FiglH 
The top diagram vanishes for SU(2) as an odd number of color charges are brought to a 
single point which vanish by color averaging. Their contribution is 
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1 1 

2 3! 

1 1 



o- a' 

^,^^P'Y.^aq2a'Y.^a'q2a'' I dfd^G\^-r)G\Q) 



4/ 



(VI.26) 



The three loop contribution stemming from three- and four-particle density is shown in 
Fig. [5l These contributions will not be quoted here. They only contribute for N^. > 2. They 
can be shown to arise from color magnetism, thus subleading in the electric cQGP under 
considerations. 




FIG. 5: Three loop contribution from three and four particle interactions 



D. Higher Loops 

In so far, the color averaging at one- two and three-loops have led to simple powers 
of Casimirs. This observation does not carry simply at higher orders as more complex 
combination of Casimirs appear. Indeed, consider the five loop contribution shown in Fig. [51 
Its contribution can be found explicitly as 

H/^' naq2a' Ea' na'q2a'^ I dfd^'G^f' - f) SU{2) 
< Ea ria Ea' Ua' (^gs'^s' - j,q3'{q2a' + q2a'') + §q2aV) (VI.27) 

X Jdfdf'G^f' -f) SU{3) 
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FIG. 6: A typical five-loop contribution. See text. 

The color averaging leads a more complex combination of Casimirs for SU(3) in comparison 
to SU(2). The cubic Casimir for SU(3) is absent for SU(2). The way the color averaging 
occurs is by noting that each vertex involves 6 lines in Fig. [61 In SU(2) the averaging 
involves only the quadratic Casimir, and the result with 3 quadratic Casimirs follow. This 
observation explains why the SU(2) color averaging in (]IV.16I) (low density expansion) can 
not be extended to SU(3) in closed form. The culpright is the occurence of the third Casimir. 



E. Result 



Now we combine the one, two and (partial) three loop results for the grand partition 
function, namely flVL23|) . NTMf and dVL26l) and use flVLT2D . to get 



^ = (E^«) + l^OhN! - i)(E^a/5C7,.)i 



+ ^47r(|-)3(iVf - ir{Y.Mf3C2ar){j:naf3Q 



2a 



I 2 



+ ^^(4vr)i(|^)5(iV,2 - l)'{Y.naiPC,^?)iT.^aPC,a 
+ ^(4vr)t(f)^iV.2-l)^(E^.(/5C^2a)f (E^«(/5^2.))^e-2- 

for both SU(2) and SU(3), with in addition for SU(3) alone 



-^|(|^)^4.r(En.)^/5^fi?i(3«:.a) 

(VI.29) 
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To proceed further, we set all classical fugacities for the three species a = 1,2,3 (quark, 
antiquark, gluon) to be the same, ria = na> = A. We also rescale all dimensions to the core 
size a to be reinstated by inspection as announced. Specifically A = cr^A, with all radial 
integrations cutoff by the core. Using flVI.281) and flVI.29|) with = |^/5C2a/o", we have for 
both SU(2) and SU(3) 



^ a a " a a 

+ ~Xl^hAn)l{N! - l)^(E^a')(E^a)^ + A^^(4vr)i(iV,^ - 1)^(E ^a')'(E ^a)^e-^-- 

-^•^ a a a a 

(VI.30) 

while in addition for SU(3), 



2 3! Mvr^ 'V 



(VI.31) 

We note that the core integrations stemming from | r — F |< a reduce to 



^ = ->^YN! - l)(E/5^2.)' / dfG^r) 

+A^^(iV.^-l)^(E(/5^2.)^)(E/5^2.) / dfG'{r)G{0) 
-X'^iN', - l)'(E(/^^2a)')(E/3^2a) / dfG'ir)G'{0) (VI.32) 
where the divergences are lumped in G(0). Our final result for SU(2) and SU(3) is then 



r3 

y AE+A^^(4vr)niVe'-l)(E^a) 



a a 



-Aii(47r)^(iVf - l)(E^a)^(l - e-^(^-)^^^E.-)^; 
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° a a 

^ a a 

X (e-"("-)^ "^^^ + 4(47r) ^ {J2 Saf^E^ (4(47r) ^ ^ J I )) 

a a 

a a 
3! O a a 

with in addition for SU(3) 

-A^~4.(|i)^i:)=(C|)^;i(3(4.)«j(i:.ji) 

-Ai~(te)'(£)Vo' - l)(Y.^.)i^(^)El{m^)i\\(Y.Qi) (VI.34) 

This our final result for up to three loops ignoring the diagrams of Fig. O The latters 
yield contributions that are of order A^ or higher. 

F. High Temperature 

The results of the high-temperature expansion can be recovered from the three loop 
results (IVI.SSp . For that we need to go back to the unintegrated contributions in r-space 
and expand in powers of (3 the exponentials. Up to order (3^ the result is 

^ =AE+A^(4^)niVf-l)(E7a)t 

*^ a "J a 

->^'\{Nl-l){Y.laf I rfrl + Atl(47r)niV.^-l)(E7a)t / dr^- 

4 jj Jr<a T Z ^ Jr<a T 

-X'\{A7:){Nl-l){Y.laf I dr + X'\A'n{Nl-l)\Y.la){Y.^a) 

)'(E)^(^) / rffV^-)^^^(S.-)^^ + (9(/5i) (VI.35) 
2d!^47r^ ^ 4 Ja<r<oo r'^ 
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where we used = 7/cr from (111.12p . For SU(2) the latter contribution is absent. This 
result agrees with fTvril) to order (3^1'^ for /ij,^^ (r) = —9 [a — r). There is a difference to order 
f3^ which can be traced back to the handling of the core potential hQ'\r) and w{r) in each of 
the two expansions. In the loop expansion we used an infinite core with e"*"^^'^ = 1 for r > a. 
In reality e""'^''^ is a function of r much like h^o\r) when r > a. A similar observation was 
made for hard sphere liquids in [9|. 



VII. FREE ENERGIES OF CQGP 

The free energy follows from the grand partition function which we have now constructed 
for both the density (virial) expansion and the loop (high) temperature expansion. The 
latter is understood for a liquid of hard spheres with hQ'\r) = —9{a — r) as in ( ]VI.35p up to 
order The former is given in (jIV.151) up to order A^. 



Specializing to SU(2) one component plasma(X]a7a 7) in (1VI.35P and carrying the 



r-integration with the hard core potential in mind as noted, we have 



-A'(N^ - 1)7V3 + X^-MN^ - 1)V + O(fji) (VII.l) 

3 2 

The density c(= A|^(^)) is 

c = A + Ac = A + V^{N^ - l^h^ + 0{(3^) (VII.2) 

The corresponding shift induced in the chemical potential, = /3/iA + /3A/ic, due to the 
interactions can be extracted from 



c = = 9p^,,+(3A^,, ^ 9p^,, ^ . _ ^ ^ A/5A/i, = A + Ac (VII.3) 

A3 A3 ^3 ^ ^ ^ 

with 



PAfx, = ^ = V^{Nl - l)7iA5 + C>(/?2) . (VII.4) 
A 

Using this shift and defining the free energy through the legendre transform 
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= In A) + (VII.5) 
we obtain for either the loop expansion or the high temperature expansion 



=-c+^(iV2-l)AM+vr(iV,2-l)A27V 

= -A - '^-^{Nl - l)7iAi - v^A^(iV,2 _ 1)^1 1^^(^2 _ 1)^1^1^ 
+ ]^T^{Nl - 1)VA' + vr(iV2 - 1)AV^ - 27ri(iV2 - l)}^^ia' + c/3/iA + 
= -A - ^-^{Nl - l)7iAi - V^xHnI - 1)7^ (Ac + 0{^')) + ^(Ac)^ 

- -A(l + ^) - ^(iV.^ - l)7iAi (1 + .(iV.^ - 1)A^(1 + 

-27ri (iV2 - l)Ai (l + ^) ^7^^^' + c(/5/iA + /?A/^,) + 



A 

3 3 



. ^^(AT^ - l)7ic5 + 7r(Ar2 - l)cVf^ - 27r5(A/'^ - 1)057^^' 
+c(3^i, + . (VII.6) 

From flVII.2p . we note that Ac ~ (3^ so that A can be substituted by the concentration 
c for terms that are of order or higher to accuracy For the free energy, the high 

temperature expansion deviates from the loop-expansion by such substitutions in higher 
order. 

The free energy in the low density expansion can be extracted from pV.15|) using a similar 



transform and substitution. For SU(2) the procedure of subtitution of A by c is detailed in 
For SU(2) D2 = B2 + fS^/^ where the B's are Bernoulli's numbers. For a one species SU(2) 
plasma e = -y/e^ = ea- Recalling that in the normalizations were carried using the 
Wigner-Seitz radius instead of the core size with 6 = cr/aws, we have e = T/6. The SU(2) 
free energy of a one- species cQGP in the low density expansion is 



Fioyj{l3,c) n n ^ ^ Q 3 2 9 3 2 /"^^sinht 

j3 — = — c — 2vvr72c2 — 37r7 c — -7r7 c / at 

V 4 Jo t 
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+ — — ( 2 sinh(— 7) + 2— cosh(— 7) + 3— sinh(— 7) + 9— cosh(-7) 



12 7 V cr a a cr^ a 6 

+ c(3^, + 0{ci) . (VII.7) 

VIII. EXCESS FREE ENERGIES FOR CQGP 

The free energies obtained above can be rewritten in terms of the plasma constant 

r = ^ (VIII. 1) 

47r Taws 

with kB = I and aws the Wigner-Seitz radius satisfying N/V{4:7!^a^rg/3) = 1. C2 is the 
quadratic Casimir, C2 = q2/{N^ — 1) and g is the strength of the coupling. Since the 
Wigner-Seitz radius is related to the density, N/V = c, it is straightforward to rewrite the 
expanded free energies above in terms of F. For instance, the first two terms in the loop 
expansion (1VII.6P read 



NT N N 3 Uti' ^ ' 

1 .^4^^i I^^C^. 

V3 ^ 3 ^ 47r aws 

= -l-^iN!-l)Tl (VIII.2) 

If we define the excess free energy Fg^ as F{T) = F{0) + Fe^(F), we get the excess free 
energies for ( jVII.61) and ( ]VII.7I) as 



F, 



loop, ex 



-^{N^ - l)Fi + h{N^ - 1)F2 - 3V36\N^ - 1)tI + 0{T^) 



NT ^ ^ ' 4 

^io^,ex ^ _y3p| _ 9^3 _ ^p3 f '^^ sinht 
NT A 16 Jo t 

+ (2 sinh(^F) + 2^ cosh(^F) + 3^ sinh(^F) + 9^ cosh(^7))i VIII.3) 

To compare these expanded results with the full molecular dynamics simulations in 
we may also write 
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mol,ex 



_4.9 r '^-2T + 3.2 ■ 4r^ - 2.2 ■ 41-3 + P--oiA+^) .yjjj 4) 
NT J+o r NT ^ ^ 

with A = 2/3(2/(f A^elA^c - 1))^)^/^. obtained from the potential energy (36) 

in [si through the following relation 



iVT J+o NT r NT 
For completeness we also quote the excess energy in the Debye-Huckel limit |8|], 



(VIII.5) 



^ = -(iV,^-l)y((Ar)itan-^(-l=)-l(ln(l + Ar)-Ar)) . (VIII.6) 




FIG. 7: Excess free energies vs Debye-Huckel (a) and molecular dynamics (b). See text. 

In Fig. [7^ we show the behavior of the excess free energies for the low density and 
loop expansions in comparison to the Debye-Huckel result. By about F ^ 1.5 the three 
expansions deviate substantially and are no longer reliable for the free energy. The deviation 



occurs earlier in the energy density as we discussed in^Ul]. In Fig. Uh we compare the three 
expansions to the molecular dynamics simulations in [3] for 0.2 < F < 1.2. We note that the 
ln(+0) contribution in flVIII.4p combines with F(+0) to zero(Infrared renormalization). The 
molecular dynamics simulations use a finite 3-box in space which at weak coupling does not 
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accomodate the Debye cloud which lowers (attractive) the free energy. Thus the discrepancy 
with the Debye-Huckel limit which is exact for a classical plasma. Simulations with larger 
3-boxes are possible but numerically time consuming. At strong coupling, the molecular 
simulations readily accomodate the short range liquid and/or crystal correlations with the 
current finite 3-boxes. They are more reliable at large F. 

IX. CONCLUSIONS 

We have explicitly analyzed the grand partition function of a cQGP defined by colored 
spheres with a hard core potential and a long range Coulomb field, borrowing on methods 
used for classical liquids 7|. The partition function was worked out explicitly for both 
SU(2) and SU(3) at low densities through a cumulant expansion, and at high temperature 
through a loop expansion. Both expansions are shown to be valid for low plasma coupling 
r = V/K ^ 1. This is the regime where the quantum QGP is expected to be dilute. Current 
interest in the quantum QGP is in the range of F 3 were the Coulomb interactions are 
much stronger and the phase liquid-like. In the follow up paper HI] we construct an equation 
of state that interpolates smoothly between the expanded results of this paper at low F and 
the strongly coupled molecular dynamics results at large F 

X. ACKNOWLEDGMENTS 

This work was supported in part by US-DOE grants DE-FG02-88ER40388 and DE-FG03- 
97ER4014. 

APPENDIX A: COLOR CHARGES 

The set of classical color charges form a Darboux's set 

with canonical Poisson brackets 

{<l>a,7Tp} = 6^p {Q'^,Q^} = r^^Q^ (A.2) 
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where is the number of colors and f^^^ is the structure constant of SU(A''c). The exphcit 
SU(2) and SU(3) representations of these sets can be found in (4, [l2|. In this Appendix, we 



quote some formulas for completeness. 
The Darboux's set for SU(2) is Q 



gi = cos0iJJ2 - vrf, = sinM - Tcf, Q 



(A.3) 



where J is used to represent the quadratic Casimir with q2 = J^a" Q^'Q"'- The phase 
space SU(2) measure used in the text is 



dQ = CRd'Kid(l)iJdJ5{j'^ — ^2) 



(A.4) 



The representation dependent constant cr is chosen so that f dQ 



The Darboux's set for SU(3) is more involved. Following 



we define the canonical set 



as (7ri,7r2, 713,01,02, 



Ji, J2) in terms of which the 8 color charges read 



with 



= cos0i7r_|_7r_ = sin 0i7r_|_7r_ = tti = tt2 

Q4 = C++71+A + C+^Tl-B = S++TT+A + S+^Tl-B 



(A.5) 



and 



71"+ = VtTs + TTi, C±± = COS ^i±(t>l + ± 03 y 

^(±01 + V^02±03: 



7r_ = V^Ts — TTi, S±± = sin 



(A.6) 



A 



B 



1 

2^V' 



J1-J2 



, vr2 Ji + 2J2 



7^2 w ^'^l + '^2 



7r2 

73 



- -^2 , 7r2 Ji + 2J2 



, 7r2 2Ji + J2 



V3 



^) (A.7) 
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Ji and J2 define the quadratic Casimir q2 and cubic Casimir ^3. Specifically, 



?2 = E = ^(^1 + J1J2 + Jl) 

?3 = E dapjQ'^Q^Q'' = Y^iJ^- •^2)(Jl + 2J2)(2Ji + J2) (A.8) 
with the help of the SU(2) symmetric tensor dap-y- we recall that da/s-y — d°'^'^ and d°'^'^ — 

To carry some of the color integrations in the text, we quote some useful identities involving 
(i's. For instance. 



CKi/?>7 

E c^^^^rf^'^^ = (A.9) 



and 



E d'^^'d'^'^'d'^^'^d^^^ = 253' 

a,/9,7,5,e,C 

E d'^^'d^^'d'^^^d'^^^ = -2g3 (A. 10) 

a,/3,7,<5,e,C 



The SU(3) phase space measure is 



\/3 

c^Q = CRd(j)id(f)2d(f)3dnidn2d7rsdJidJ2—JiJ2{Ji + J2) 

4o 



X5(|(J2 + Jl J2 + J|) - ?2)5(^(Jl - J2)(Jl + 2J2)(2Ji + J2) - ?3) (A.ll) 

with again cr set by the normalization of the color space volume to 1. 



APPENDIX B: COLOR INTEGRATIONS 

In this Appendix we explicit some of the color integrations carried in the text. For SU(2) 
all color charge integrations can be done analytically. The representation dependent cr is set 
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TABLE I: Non-zero constants of d"^'^ for SU(3) 

^118 ^228 ^338 ^448 ^558 ^668 ^778 ^888 

J_ J_ J_ 1 1 1 1 ]_ 

2^3 2^3 2^3 2V3 Vs 

^146 ^157 ^247 ^256 ^344 ^355 ^366 ^377 

1 1 _1 1 1 1 _1 _1 

2222 2 2 22 

to give the normalization of the colored space volume to be 1. For instance c/j = l/27ry/g^[l/ 
by using / dQ = 1. The SU(3) integrations cannot be done analytically given the cubic nature 
of the constraint in the color measure. The quadratic and cubic Casimirs are fixed by 

J dQq2 = J dQ^iJ^ + J1J2 + J^) = q2 

J dQqs = J dQ^iJi - J2)(Ji + 2J2)(2Ji + J2) = qs (B.l) 

One Color: 



rfQQ" = . (B.2) 



Two Colors: 



J dQQ'^Q^ = ^2^^ = ^ylry^"^ (B.3) 

with C2 = Nc for gluons and C2 = 1/2 for quarks. For SU(3), the cubic Casimir q-s = for 
gluons and gs = (A^^^ - l)(A^c ~ 4)/4A^c for quarks. It is zero for SU(2). 
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Three Colors: 



E d^p^Q'^Q^Q'^ = 53 = i ^ = j dQq, = E d^P-y J dgQ^^Q^Q^ (B.4) 

i,/3,7 c*i/3i7 Cf,/3,7 



where we used ( lA.Sp and (IB.ip and the identity 



E rf-^^rf.^, = — (iV,2 - l)(iVf - 4) = 4g3 



(B.5) 



Thus 



(B.6) 



Four Colors: 



/ dQQ'^Q^Q^Q^ = A{Y,d"^'^d^^'^ + Y,d"^'^d'^^'^ + Y.d"^'^d^^'^) + B{5''^5^^ + 5"'<5'^^ + 5"^5^'<) 

n n n 

(B.7) 

For SU(2) A = and B = gl/lS. For SU(3) this integral is more involved. However, for 
certain arrangements of charges, say Q^,Q^,Q^,Q^ the integrations can be done and thus 
the constants A, B fixed. For that we need to undo integrals involving tt's since = tt"^. 
The TT integrations are bounded 0, [l^ . To carry them, it is best to change variables 



so that 



dTT2d'K^ = —^dxdy (B.9) 



It is also useful to define 
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so that 



Ki = -(2Ji + J2), K2 = -(Ji + 2J2) 



(B.IO) 



K2-Ki<x<Ki, Ki-K2<y<K2 
With these definitions, we can esily unwind some vr integrations. For instance 



(B.ll) 



and 



diXid-K 2(171 z = — — / dx dy drci = —-J1J2U1 + J2) 

2 JK2-K1 Jk^-K2 J-U^+y) 4 



d'Ki'K2'd'K2d'K-i 



^ r dx 

2 JK2-K1 



dyC^ix - y)f ^ dTXi 



l^JMJi + J2)^(Ji' + J1J2 + Jl) . 



This latter integral is directly related to 



(B.12) 



(B.13) 



j dQQ^Q^ = cr J d(f)id(f)2d(f)3 j dT[i'K2^dT[2dT[-i j dJidJ2^JiJ2{Ji + J2) 

x5{\{Jl + J1J2 + Jl) - g2)5(-^(Ji - J2){Ji + 2J2)(2Ji + J2) - qs) 
6 io 



= / ^^^^(^1 + "^1^2 + Jl) = J dQ^q2 = j;^^q2 
after using flA.lip and fIB.ll) . The result agrees with flB.3p . Similarly, 



(B.14) 



dTliTX2^dTX2dll3 



2 JK2-K1 JK1-K2 



dy{^{x-y) 



ii^+y) 



diTi 



^^^1^2(^1 + J2)-{Ji - J2){Ji + 2J2)(2Ji + J2) (B.15) 



is related to 



f f f 

dQQ^Q^Q^ = crJ d(f)id(f)2d(f)s J d7ii7r2^dn2d7r3 J dJidJ2 — JiJ2{Ji + J-. 
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x5(^(j2 + J1J2 + J|) - q2)6i^iJ^ - J2){Ji + 2J2)(2Ji + J2) - Qz) 



40 

f ,^ 113 1 



J rfQ^(Ji - J2)(Ji + 2J2)(2Ji + J2) 



--loi^^^^^-iTIT^^^^i' (B.16) 

after using = {N^ — 1){N^ — 4:)/4:Nc = 10/3 for SU(3). This result is in agreement with 
( IB. 61) . A simihrar reasoning, yields 

A + 3B = ^ql (B.17) 

for the A, B constants in fIB.Tp . 



Five Colors: 



J 560 ^ 

^^/37^a5e ^ ^/35^a7e ^ d'^^S + 5^^ d""^' + + d^'d^^^A (B.18) 



with 



^ :g2g3 = Wq2^q-i = l;C2 (B.19) 



560" " 78" 10" 7 



Six Colors: 
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(B.20) 

with the A, B constants tabulated below. 

TABLE II: Constants A and B in the integration of six color charges 



A B 



SU{2) 
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